Minimum information and guidelines for reporting a multiplexed assay of variant effect.

Multiplexed assays of variant effect (MAVEs) have emerged as a powerful approach for interrogating thousands of genetic variants in a single experiment. The flexibility and widespread adoption of these techniques across diverse disciplines have led to a heterogeneous mix of data formats and descriptions, which complicates the downstream use of the resulting datasets. To address these issues and promote reproducibility and reuse of MAVE data, we define a set of minimum information standards for MAVE data and metadata and outline a controlled vocabulary aligned with established biomedical ontologies for describing these experimental designs.

The electrophysiologic effects of KCNQ1 extend beyond expression of IKs: evidence from genetic and pharmacologic block.

AIMS: While variants in KCNQ1 are the commonest cause of the congenital long QT syndrome, we and others find only a small IKs in cardiomyocytes from human induced pluripotent stem cells (iPSC-CMs) or human ventricular myocytes. METHODS AND RESULTS: We studied population control iPSC-CMs and iPSC-CMs from a patient with Jervell and Lange-Nielsen (JLN) syndrome due to compound heterozygous loss of function KCNQ1 variants. We compared the effects of pharmacologic IKs block to those of genetic KCNQ1 ablation, using JLN cells, cells homozygous for the KCNQ1 loss of function allele G643S, or siRNAs reducing KCNQ1 expression. We also studied the effects of two blockers of IKr, the other major cardiac repolarizing current, in the setting of pharmacologic or genetic ablation of KCNQ1: moxifloxacin, associated with a very low risk of drug-induced long QT, and dofetilide, a high-risk drug.In control cells, a small IKs was readily recorded but pharmacologic IKs block produced no change in action potential duration at 90% repolarization (APD90). By contrast, in cells with genetic ablation of KCNQ1 (JLN), baseline APD90 was markedly prolonged compared with control cells (469 ± 20 vs. 310 ± 16 ms). JLN cells displayed increased sensitivity to acute IKr block: the concentration (µM) of moxifloxacin required to prolong APD90 100 msec was 237.4 (median, IQR 100.6-391.6, n = 7) in population cells versus 23.7 (17.3-28.7, n = 11) in JLN cells. In control cells, chronic moxifloxacin exposure (300µM) mildly prolonged APD90 (10%) and increased IKs, while chronic exposure to dofetilide (5 nM) produced greater prolongation (67%) and no increase in IKs. However, in the siRNA-treated cells, moxifloxacin did not increase IKs, and markedly prolonged APD90. CONCLUSION: Our data strongly suggest that KCNQ1 expression modulates baseline cardiac repolarization, and the response to IKr block, through mechanisms beyond simply generating IKs. TRANSLATIONAL PERSPECTIVE: Mutations in KCNQ1 - whose expression generates IKs - are the major cause of long QT syndrome. We report here that while pharmacologic IKs block in human cardiomyocytes generates minimal change in repolarization, suppressing KCNQ1 expression markedly increases both baseline repolarization duration and sensitivity to some (but not all) specific IKr blockers. Thus, beyond simply generating IKs, KCNQ1 subserves critical additional role(s) in repolarization control at baseline and in response to IKr block. Our findings imply that assessment of arrhythmic risk in individual patients and by drugs requires a framework that extends beyond a simple one gene-one ion current paradigm.

Multifocal Ectopic Purkinje Premature Contractions due to neutralization of an SCN5A negative charge: structural insights into the gating pore hypothesis.

Background: We identified a novel SCN5A variant, E171Q, in a neonate with very frequent ectopy and reduced ejection fraction which normalized after arrhythmia suppression by flecainide. This clinical picture is consistent with multifocal ectopic Purkinje-related premature contractions (MEPPC). Most previous reports of MEPPC have implicated SCN5A variants such as R222Q that neutralize positive charges in the S4 voltage sensor helix of the channel protein NaV1.5 and generate a gating pore current. Methods and Results: E171 is a highly conserved negatively-charged residue located in the S2 transmembrane helix of NaV1.5 domain I. E171 is a key component of the Gating Charge Transfer Center, a region thought to be critical for normal movement of the S4 voltage sensor helix. We used heterologous expression, CRISPR-edited induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs), and molecular dynamics simulations to demonstrate that E171Q generates a gating pore current, which was suppressed by a low concentration of flecainide (IC50 = 0.71±0.07 µM). R222Q shifts voltage dependence of activation and inactivation in a negative direction but we observed positive shifts with E171Q. E171Q iPSC-CMs demonstrated abnormal spontaneous activity and prolonged action potentials. Molecular dynamics simulations revealed that both R222Q and E171Q proteins generate a water-filled permeation pathway that underlies generation of the gating pore current. Conclusion: Previously identified MEPPC-associated variants that create gating pore currents are located in positively-charged residues in the S4 voltage sensor and generate negative shifts in the voltage dependence of activation and inactivation. We demonstrate that neutralizing a negatively charged S2 helix residue in the Gating Charge Transfer Center generates positive shifts but also create a gating pore pathway. These findings implicate the gating pore pathway as the primary functional and structural determinant of MEPPC and widen the spectrum of variants that are associated with gating pore-related disease in voltage-gated ion channels.

Prognostic Value of Multiplexed Assays of Variant Effect and Automated Patch-clamping for KCNH2-LQTS Risk Stratification.

Background: Long QT syndrome (LQTS) is a lethal arrhythmia condition, frequently caused by rare loss-of-function variants in the cardiac potassium channel encoded by KCNH2. Variant-based risk stratification is complicated by heterogenous clinical data, incomplete penetrance, and low-throughput functional data. Objective: To test the utility of variant-specific features, including high-throughput functional data, to predict cardiac events among KCNH2 variant heterozygotes. Methods: We quantified cell-surface trafficking of 18,323 variants in KCNH2 and recorded potassium current densities for 506 KCNH2 variants. Next, we deeply phenotyped 1150 KCNH2 missense variant patients, including ECG features, cardiac event history (528 total cardiac events), and mortality. We then assessed variant functional, in silico, structural, and LQTS penetrance data to stratify event-free survival for cardiac events in the study cohort. Results: Variant-specific current density (HR 0.28 [0.13-0.60]) and estimates of LQTS penetrance incorporating MAVE data (HR 3.16 [1.59-6.27]) were independently predictive of severe cardiac events when controlling for patient-specific features. Risk prediction models incorporating these data significantly improved prediction of 20 year cardiac events (AUC 0.79 [0.75-0.82]) over patient-only covariates (QTc and sex) (AUC 0.73 [0.70-0.77]). Conclusion: We show that high-throughput functional data, and other variant-specific features, meaningfully contribute to both diagnosis and prognosis of a clinically actionable monogenic disease.

ParSE-seq: A Calibrated Multiplexed Assay to Facilitate the Clinical Classification of Putative Splice-altering Variants.

Background: Interpreting the clinical significance of putative splice-altering variants outside 2-base pair canonical splice sites remains difficult without functional studies. Methods: We developed Parallel Splice Effect Sequencing (ParSE-seq), a multiplexed minigene-based assay, to test variant effects on RNA splicing quantified by high-throughput sequencing. We studied variants in SCN5A, an arrhythmia-associated gene which encodes the major cardiac voltage-gated sodium channel. We used the computational tool SpliceAI to prioritize exonic and intronic candidate splice variants, and ClinVar to select benign and pathogenic control variants. We generated a pool of 284 barcoded minigene plasmids, transfected them into Human Embryonic Kidney (HEK293) cells and induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs), sequenced the resulting pools of splicing products, and calibrated the assay to the American College of Medical Genetics and Genomics scheme. Variants were interpreted using the calibrated functional data, and experimental data were compared to SpliceAI predictions. We further studied some splice-altering missense variants by cDNA-based automated patch clamping (APC) in HEK cells and assessed splicing and sodium channel function in CRISPR-edited iPSC-CMs. Results: ParSE-seq revealed the splicing effect of 224 SCN5A variants in iPSC-CMs and 244 variants in HEK293 cells. The scores between the cell types were highly correlated (R2=0.84). In iPSCs, the assay had concordant scores for 21/22 benign/likely benign and 24/25 pathogenic/likely pathogenic control variants from ClinVar. 43/112 exonic variants and 35/70 intronic variants with determinate scores disrupted splicing. 11 of 42 variants of uncertain significance were reclassified, and 29 of 34 variants with conflicting interpretations were reclassified using the functional data. SpliceAI computational predictions correlated well with experimental data (AUC = 0.96). We identified 20 unique SCN5A missense variants that disrupted splicing, and 2 clinically observed splice-altering missense variants of uncertain significance had normal function when tested with the cDNA-based APC assay. A splice-altering intronic variant detected by ParSE-seq, c.1891-5C>G, also disrupted splicing and sodium current when introduced into iPSC-CMs at the endogenous locus by CRISPR editing. Conclusions: ParSE-seq is a calibrated, multiplexed, high-throughput assay to facilitate the classification of candidate splice-altering variants.

Minimum information and guidelines for reporting a Multiplexed Assay of Variant Effect.

Multiplexed Assays of Variant Effect (MAVEs) have emerged as a powerful approach for interrogating thousands of genetic variants in a single experiment. The flexibility and widespread adoption of these techniques across diverse disciplines has led to a heterogeneous mix of data formats and descriptions, which complicates the downstream use of the resulting datasets. To address these issues and promote reproducibility and reuse of MAVE data, we define a set of minimum information standards for MAVE data and metadata and outline a controlled vocabulary aligned with established biomedical ontologies for describing these experimental designs.

High-throughput functional mapping of variants in an arrhythmia gene, KCNE1, reveals novel biology.

Background: KCNE1 encodes a 129-residue cardiac potassium channel (IKs) subunit. KCNE1 variants are associated with long QT syndrome and atrial fibrillation. However, most variants have insufficient evidence of clinical consequences and thus limited clinical utility. Results: Here, we demonstrate the power of variant effect mapping, which couples saturation mutagenesis with high-throughput sequencing, to ascertain the function of thousands of protein coding KCNE1 variants. We comprehensively assayed KCNE1 variant cell surface expression (2,554/2,709 possible single amino acid variants) and function (2,539 variants). We identified 470 loss-of-surface expression and 588 loss-of-function variants. Out of the 588 loss-of-function variants, only 155 had low cell surface expression. The latter half of the protein is dispensable for protein trafficking but essential for channel function. 22 of the 30 KCNE1 residues (73%) highly intolerant of variation were in predicted close contact with binding partners KCNQ1 or calmodulin. Our data were highly concordant with gold standard electrophysiological data (ρ = -0.65), population and patient cohorts (32/38 concordant variants), and computational metrics (ρ = -0.55). Our data provide moderate-strength evidence for the ACMG/AMP functional criteria for benign and pathogenic variants. Conclusions: Comprehensive variant effect maps of KCNE1 can both provide insight into IKs channel biology and help reclassify variants of uncertain significance.

Multicenter clinical and functional evidence reclassifies a recurrent noncanonical filamin C splice-altering variant.

BACKGROUND: Truncating variants in filamin C (FLNC) can cause arrhythmogenic cardiomyopathy (ACM) through haploinsufficiency. Noncanonical splice-altering variants may contribute to this phenotype. OBJECTIVE: The purpose of this study was to investigate the clinical and functional consequences of a recurrent FLNC intronic variant of uncertain significance (VUS), c.970-4A>G. METHODS: Clinical data in 9 variant heterozygotes from 4 kindreds were obtained from 5 tertiary health care centers. We used in silico predictors and functional studies with peripheral blood and patient-specific induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs). Isolated RNA was studied by reverse transcription polymerase chain reaction. iPSC-CMs were further characterized at baseline and after nonsense-mediated decay (NMD) inhibition, using quantitative polymerase chain reaction (qPCR), RNA-sequencing, and cellular electrophysiology. American College of Medical Genetics and Genomics (ACMG) criteria were used to adjudicate variant pathogenicity. RESULTS: Variant heterozygotes displayed a spectrum of disease phenotypes, spanning from mild ventricular dysfunction with palpitations to severe ventricular arrhythmias requiring device shocks or progressive cardiomyopathy requiring heart transplantation. Consistent with in silico predictors, the c.970-4A>G FLNC variant activated a cryptic splice acceptor site, introducing a 3-bp insertion containing a premature termination codon. NMD inhibition upregulated aberrantly spliced transcripts by qPCR and RNA-sequencing. Patch clamp studies revealed irregular spontaneous action potentials, increased action potential duration, and increased sodium late current in proband-derived iPSC-CMs. These findings fulfilled multiple ACMG criteria for pathogenicity. CONCLUSION: Clinical, in silico, and functional evidence support the prediction that the intronic c.970-4A>G VUS disrupts splicing and drives ACM, enabling reclassification from VUS to pathogenic.

Transcriptional Dysregulation Underlies Both Monogenic Arrhythmia Syndrome and Common Modifiers of Cardiac Repolarization.

BACKGROUND: Brugada syndrome (BrS) is an inherited arrhythmia syndrome caused by loss-of-function variants in the cardiac sodium channel gene SCN5A (sodium voltage-gated channel alpha subunit 5) in ≈20% of subjects. We identified a family with 4 individuals diagnosed with BrS harboring the rare G145R missense variant in the cardiac transcription factor TBX5 (T-box transcription factor 5) and no SCN5A variant. METHODS: We generated induced pluripotent stem cells (iPSCs) from 2 members of a family carrying TBX5-G145R and diagnosed with Brugada syndrome. After differentiation to iPSC-derived cardiomyocytes (iPSC-CMs), electrophysiologic characteristics were assessed by voltage- and current-clamp experiments (n=9 to 21 cells per group) and transcriptional differences by RNA sequencing (n=3 samples per group), and compared with iPSC-CMs in which G145R was corrected by CRISPR/Cas9 approaches. The role of platelet-derived growth factor (PDGF)/phosphoinositide 3-kinase (PI3K) pathway was elucidated by small molecule perturbation. The rate-corrected QT (QTc) interval association with serum PDGF was tested in the Framingham Heart Study cohort (n=1893 individuals). RESULTS: TBX5-G145R reduced transcriptional activity and caused multiple electrophysiologic abnormalities, including decreased peak and enhanced "late" cardiac sodium current (INa), which were entirely corrected by editing G145R to wild-type. Transcriptional profiling and functional assays in genome-unedited and -edited iPSC-CMs showed direct SCN5A down-regulation caused decreased peak INa, and that reduced PDGF receptor (PDGFRA [platelet-derived growth factor receptor α]) expression and blunted signal transduction to PI3K was implicated in enhanced late INa. Tbx5 regulation of the PDGF axis increased arrhythmia risk due to disruption of PDGF signaling and was conserved in murine model systems. PDGF receptor blockade markedly prolonged normal iPSC-CM action potentials and plasma levels of PDGF in the Framingham Heart Study were inversely correlated with the QTc interval (P<0.001). CONCLUSIONS: These results not only establish decreased SCN5A transcription by the TBX5 variant as a cause of BrS, but also reveal a new general transcriptional mechanism of arrhythmogenesis of enhanced late sodium current caused by reduced PDGF receptor-mediated PI3K signaling.

Multi-site validation of a functional assay to adjudicate SCN5A Brugada Syndrome-associated variants.

Brugada Syndrome (BrS) is an inheritable arrhythmia condition that is associated with rare, loss-of-function variants in the cardiac sodium channel gene, SCN5A. Interpreting the pathogenicity of SCN5A missense variants is challenging and ~79% of SCN5A missense variants in ClinVar are currently classified as Variants of Uncertain Significance (VUS). An in vitro SCN5A-BrS automated patch clamp assay was generated for high-throughput functional studies of NaV1.5. The assay was independently studied at two separate research sites - Vanderbilt University Medical Center and Victor Chang Cardiac Research Institute - revealing strong correlations, including peak INa density (R2=0.86). The assay was calibrated according to ClinGen Sequence Variant Interpretation recommendations using high-confidence variant controls (n=49). Normal and abnormal ranges of function were established based on the distribution of benign variant assay results. The assay accurately distinguished benign controls (24/25) from pathogenic controls (23/24). Odds of Pathogenicity values derived from the experimental results yielded 0.042 for normal function (BS3 criterion) and 24.0 for abnormal function (PS3 criterion), resulting in up to strong evidence for both ACMG criteria. The calibrated assay was then used to study SCN5A VUS observed in four families with BrS and other arrhythmia phenotypes associated with SCN5A loss-of-function. The assay revealed loss-of-function for three of four variants, enabling reclassification to likely pathogenic. This validated APC assay provides clinical-grade functional evidence for the reclassification of current VUS and will aid future SCN5A-BrS variant classification.

Genetics of congenital arrhythmia syndromes: the challenge of variant interpretation.

Congenital arrhythmia syndromes are rare genetic disorders that can cause a high risk of sudden cardiac death. Expert panels have affirmed 15 genes that are linked to congenital arrhythmias. These genes mostly encode cardiac ion channel proteins or associated regulatory proteins that generate the cardiac action potential. Common genetic variation modulates the risk of rare variants and partially explains the incomplete penetrance of these disorders. As genetic testing becomes more prevalent, a major challenge is that most detected variants are annotated as variants of uncertain significance. This review will highlight emerging methods that are refining our understanding of arrhythmia genetics, including phenotype risk scores, large cohorts, in vitro functional assays, structural models, and computational predictions.

Functional Assays Reclassify Suspected Splice-Altering Variants of Uncertain Significance in Mendelian Channelopathies.

BACKGROUND: Rare protein-altering variants in SCN5A, KCNQ1, and KCNH2 are major causes of Brugada syndrome and the congenital long QT syndrome. While splice-altering variants lying outside 2-bp canonical splice sites can cause these diseases, their role remains poorly described. We implemented 2 functional assays to assess 12 recently reported putative splice-altering variants of uncertain significance and 1 likely pathogenic variant without functional data observed in Brugada syndrome and long QT syndrome probands. METHODS: We deployed minigene assays to assess the splicing consequences of 10 variants. Three variants incompatible with the minigene approach were introduced into control induced pluripotent stem cells by CRISPR genome editing. We differentiated cells into induced pluripotent stem cell-derived cardiomyocytes and studied splicing outcomes by reverse transcription-polymerase chain reaction. We used the American College of Medical Genetics and Genomics functional assay criteria (PS3/BS3) to reclassify variants. RESULTS: We identified aberrant splicing, with presumed disruption of protein sequence, in 8/10 variants studied using the minigene assay and 1/3 studied in induced pluripotent stem cell-derived cardiomyocytes. We reclassified 8 variants of uncertain significance to likely pathogenic, 1 variant of uncertain significance to likely benign, and 1 likely pathogenic variant to pathogenic. CONCLUSIONS: Functional assays reclassified splice-altering variants outside canonical splice sites in Brugada Syndrome- and long QT syndrome-associated genes.

Mortality Among Patients With Early-Onset Atrial Fibrillation and Rare Variants in Cardiomyopathy and Arrhythmia Genes.

Importance: Patients with early-onset atrial fibrillation (AF) are enriched for rare variants in cardiomyopathy and arrhythmia genes. The clinical significance of these rare variants in patients with early-onset AF is unknown. Objective: To assess the association between rare variants in cardiomyopathy and arrhythmia genes detected in patients with early-onset AF and time to death. Design, Setting, and Participants: This prospective cohort study included participants with AF diagnosed before 66 years of age who underwent whole-genome sequencing through the National Heart, Lung and Blood Institute's Trans-Omics for Precision Medicine program. Participants were enrolled from November 23, 1999, to June 2, 2015. Data were analyzed from February 26 to September 19, 2021. Exposures: Rare variants identified in a panel of 145 genes that are included in cardiomyopathy and arrhythmia panels used by commercial clinical genetic testing laboratories. Main Outcomes and Measures: The primary study outcome was time to death and was adjudicated from medical records and the National Death Index. Multivariable Cox proportional hazards regression was used to evaluate the association of disease-associated variants with risk of death after adjustment for age at AF diagnosis, sex, race, body mass index, left ventricular ejection fraction, and an interaction term of age at AF diagnosis and disease-associated variant status. Results: Among 1293 participants (934 [72%] male; median age at enrollment, 56.0 years; IQR, 48.0-61.0 years), disease-associated (pathogenic or likely pathogenic) rare variants were found in 131 (10%). During a median follow-up of 9.9 years (IQR, 6.9-13.2 years), 219 participants (17%) died. In univariable analysis, disease-associated variants were associated with an increased risk of mortality (hazard ratio, [HR], 1.5; 95% CI, 1.0-2.1; P = .05); the association remained significant in multivariable modeling when adjusted for age at AF diagnosis, sex, race, body mass index, left ventricular ejection fraction, and an interaction term between disease-associated variant status and age at AF diagnosis. The interaction demonstrated that disease-associated variants were associated with a significantly higher risk of mortality compared with no disease-associated variant when AF was diagnosed at a younger age (P = .008 for interaction). Higher body mass index (per IQR: HR, 1.4; 95% CI, 1.2-1.6; P < .001) and lower left ventricular ejection fraction (per IQR: HR, 0.8; 95% CI, 0.7-0.8; P < .001) were associated with higher mortality risk. There were 73 cardiomyopathy-related deaths, 40 sudden deaths, and 10 stroke-related deaths. Mortality among patients with the most prevalent genes with disease-associated variants was 26% (10 of 38 patients) for TTN, 33% (6 of 18) for MYH7, 22% (2 of 9) for LMNA, 0% (0 of 10) for MYH6, and 0% (0 of 8) for KCNQ1. Conclusions and Relevance: The findings suggest that rare variants in cardiomyopathy and arrhythmia genes may be associated with increased risk of mortality among patients with early-onset AF, especially those diagnosed at a younger age. Genetic testing may provide important prognostic information for patients with early-onset AF.

Dominant negative effects of SCN5A missense variants.

PURPOSE: Up to 30% of patients with Brugada syndrome (BrS) carry loss-of-function (LoF) variants in the cardiac sodium channel gene SCN5A encoding for the protein NaV1.5. Recent studies suggested that NaV1.5 can dimerize, and some variants exert dominant negative effects. In this study, we sought to explore the generality of missense variant NaV1.5 dominant negative effects and their clinical severity. METHODS: We identified 35 LoF variants (<10% of wild type [WT] peak current) and 15 partial LoF variants (10%-50% of WT peak current) that we assessed for dominant negative effects. SCN5A variants were studied in HEK293T cells, alone or in heterozygous coexpression with WT SCN5A using automated patch clamp. To assess the clinical risk, we compared the prevalence of dominant negative vs putative haploinsufficient (frameshift, splice, or nonsense) variants in a BrS consortium and the Genome Aggregation Database population database. RESULTS: In heterozygous expression with WT, 32 of 35 LoF and 6 of 15 partial LoF variants showed reduction to <75% of WT-alone peak current, showing a dominant negative effect. Individuals with dominant negative LoF variants had an elevated disease burden compared with the individuals with putative haploinsufficient variants (2.7-fold enrichment in BrS cases, P = .019). CONCLUSION: Most SCN5A missense LoF variants exert a dominant negative effect. This class of variant confers an especially high burden of BrS.

Veratridine Can Bind to a Site at the Mouth of the Channel Pore at Human Cardiac Sodium Channel NaV1.5.

The cardiac sodium ion channel (NaV1.5) is a protein with four domains (DI-DIV), each with six transmembrane segments. Its opening and subsequent inactivation results in the brief rapid influx of Na+ ions resulting in the depolarization of cardiomyocytes. The neurotoxin veratridine (VTD) inhibits NaV1.5 inactivation resulting in longer channel opening times, and potentially fatal action potential prolongation. VTD is predicted to bind at the channel pore, but alternative binding sites have not been ruled out. To determine the binding site of VTD on NaV1.5, we perform docking calculations and high-throughput electrophysiology experiments in the present study. The docking calculations identified two distinct binding regions. The first site was in the pore, close to the binding site of NaV1.4 and NaV1.5 blocking drugs in experimental structures. The second site was at the "mouth" of the pore at the cytosolic side, partly solvent-exposed. Mutations at this site (L409, E417, and I1466) had large effects on VTD binding, while residues deeper in the pore had no effect, consistent with VTD binding at the mouth site. Overall, our results suggest a VTD binding site close to the cytoplasmic mouth of the channel pore. Binding at this alternative site might indicate an allosteric inactivation mechanism for VTD at NaV1.5.

Common Ancestry-Specific Ion Channel Variants Predispose to Drug-Induced Arrhythmias.

BACKGROUND: Multiple reports associate the cardiac sodium channel gene (SCN5A) variants S1103Y and R1193Q with type 3 congenital long QT syndrome and drug-induced long QT syndrome. These variants are too common in ancestral populations to be highly arrhythmogenic at baseline, however: S1103Y allele frequency is 8.1% in African Americans and R1193Q 6.1% in East Asians. R1193Q is known to increase late sodium current (INa-L) in cardiomyocytes derived from induced pluripotent stem cells but the role of these variants in modulating repolarization remains poorly understood. METHODS: We determined the effect of S1103Y on QT intervals among African-American participants in a large electronic health record. Using cardiomyocytes derived from induced pluripotent stem cells carrying naturally occurring or genome-edited variants, we studied action potential durations (APDs) at baseline and after challenge with the repolarizing potassium current (IKr) blocker dofetilide and INa-L and IKr at baseline. RESULTS: In 1479 African-American participants with no confounding medications or diagnoses of heart disease, QT intervals in S1103Y carriers was no different from that in noncarriers. Baseline APD was no different in cells expressing the Y allele (SY, YY cells) compared with isogenic cells with the reference allele (SS cells). However, INa-L was increased in SY and YY cells and the INa-L blocker GS967 shortened APD in SY/YY but not SS cells (P<0.001). IKr was increased almost 2-fold in SY/YY cells compared with SS cells (tail current: 0.66±0.1 versus 1.2±0.1 pA/pF; P<0.001). Dofetilide challenge prolonged APD at much lower concentrations in SY (4.1 nmol/L [interquartile range, 1.5-9.3]; n=11) and YY (4.2 nmol/L [1.7-5.0]; n=5) than in SS cells (249 nmol/L [22.3-2905]; n=14; P<0.001 and P<0.01, respectively) and elicited afterdepolarizations in 8/16 SY/YY cells but only in 1/14 SS cells. R1193Q cells similarly displayed no difference in baseline APD but increased IKr and increased dofetilide sensitivity. CONCLUSIONS: These common ancestry-specific variants do not affect baseline repolarization, despite generating increased INa-L. We propose that increased IKr serves to maintain normal repolarization but increases the risk of manifest QT prolongation with IKr block in variant carriers. Our findings emphasize the need for inclusion of diverse populations in the study of adverse drug reactions.

Arrhythmia Variant Associations and Reclassifications in the eMERGE-III Sequencing Study.

BACKGROUND: Sequencing Mendelian arrhythmia genes in individuals without an indication for arrhythmia genetic testing can identify carriers of pathogenic or likely pathogenic (P/LP) variants. However, the extent to which these variants are associated with clinically meaningful phenotypes before or after return of variant results is unclear. In addition, the majority of discovered variants are currently classified as variants of uncertain significance, limiting clinical actionability. METHODS: The eMERGE-III study (Electronic Medical Records and Genomics Phase III) is a multicenter prospective cohort that included 21 846 participants without previous indication for cardiac genetic testing. Participants were sequenced for 109 Mendelian disease genes, including 10 linked to arrhythmia syndromes. Variant carriers were assessed with electronic health record-derived phenotypes and follow-up clinical examination. Selected variants of uncertain significance (n=50) were characterized in vitro with automated electrophysiology experiments in HEK293 cells. RESULTS: As previously reported, 3.0% of participants had P/LP variants in the 109 genes. Herein, we report 120 participants (0.6%) with P/LP arrhythmia variants. Compared with noncarriers, arrhythmia P/LP carriers had a significantly higher burden of arrhythmia phenotypes in their electronic health records. Fifty-four participants had variant results returned. Nineteen of these 54 participants had inherited arrhythmia syndrome diagnoses (primarily long-QT syndrome), and 12 of these 19 diagnoses were made only after variant results were returned (0.05%). After in vitro functional evaluation of 50 variants of uncertain significance, we reclassified 11 variants: 3 to likely benign and 8 to P/LP. CONCLUSIONS: Genome sequencing in a large population without indication for arrhythmia genetic testing identified phenotype-positive carriers of variants in congenital arrhythmia syndrome disease genes. As the genomes of large numbers of people are sequenced, the disease risk from rare variants in arrhythmia genes can be assessed by integrating genomic screening, electronic health record phenotypes, and in vitro functional studies. REGISTRATION: URL: https://www. CLINICALTRIALS: gov; Unique identifier; NCT03394859.

Early-Onset Atrial Fibrillation and the Prevalence of Rare Variants in Cardiomyopathy and Arrhythmia Genes.

Importance: Early-onset atrial fibrillation (AF) can be the initial manifestation of a more serious underlying inherited cardiomyopathy or arrhythmia syndrome. Objective: To examine the results of genetic testing for early-onset AF. Design, Setting, and Participants: This prospective, observational cohort study enrolled participants from an academic medical center who had AF diagnosed before 66 years of age and underwent whole genome sequencing through the National Heart, Lung, and Blood Institute's Trans-Omics for Precision Medicine program. Participants were enrolled from November 23, 1999, to June 2, 2015. Data analysis was performed from October 24, 2020, to March 11, 2021. Exposures: Rare variants identified in a panel of 145 genes that are included on cardiomyopathy and arrhythmia panels used by commercial clinical genetic testing laboratories. Main Outcomes and Measures: Sequencing data were analyzed using an automated process followed by manual review by a panel of independent, blinded reviewers. The primary outcome was classification of rare variants using American College of Medical Genetics and Genomics criteria: benign, likely benign, variant of undetermined significance, likely pathogenic, or pathogenic. Disease-associated variants were defined as pathogenic/likely pathogenic variants in genes associated with autosomal dominant or X-linked dominant disorders. Results: Among 1293 participants (934 [72.2%] male; median [interquartile range] age at enrollment, 56 [48-61] years; median [interquartile range] age at AF diagnosis, 50 [41-56] years), genetic testing identified 131 participants (10.1%) with a disease-associated variant, 812 (62.8%) with a variant of undetermined significance, 92 (7.1%) as heterozygous carriers for an autosomal recessive disorder, and 258 (20.0%) with no suspicious variant. The likelihood of a disease-associated variant was highest in participants with AF diagnosed before the age of 30 years (20 of 119 [16.8%; 95% CI, 10.0%-23.6%]) and lowest after the age of 60 years (8 of 112 [7.1%; 95% CI, 2.4%-11.9%]). Disease-associated variants were more often associated with inherited cardiomyopathy syndromes compared with inherited arrhythmias. The most common genes were TTN (n = 38), MYH7 (n = 18), MYH6 (n = 10), LMNA (n = 9), and KCNQ1 (n = 8). Conclusions and Relevance: In this cohort study, genetic testing identified a disease-associated variant in 10% of patients with early-onset AF (the percentage was higher if diagnosed before the age of 30 years and lower if diagnosed after the age of 60 years). Most pathogenic/likely pathogenic variants are in genes associated with cardiomyopathy. These results support the use of genetic testing in early-onset AF.

Estimating the Posttest Probability of Long QT Syndrome Diagnosis for Rare KCNH2 Variants.

BACKGROUND: The proliferation of genetic profiling has revealed many associations between genetic variations and disease. However, large-scale phenotyping efforts in largely healthy populations, coupled with DNA sequencing, suggest variants currently annotated as pathogenic are more common in healthy populations than previously thought. In addition, novel and rare variants are frequently observed in genes associated with disease both in healthy individuals and those under suspicion of disease. This raises the question of whether these variants can be useful predictors of disease. To answer this question, we assessed the degree to which the presence of a variant in the cardiac potassium channel gene KCNH2 was diagnostically predictive for the autosomal dominant long QT syndrome. METHODS: We estimated the probability of a long QT diagnosis given the presence of each KCNH2 variant using Bayesian methods that incorporated variant features such as changes in variant function, protein structure, and in silico predictions. We call this estimate the posttest probability of disease. Our method was applied to over 4000 individuals heterozygous for 871 missense or in-frame insertion/deletion variants in KCNH2 and validated against a separate international cohort of 933 individuals heterozygous for 266 missense or in-frame insertion/deletion variants. RESULTS: Our method was well-calibrated for the observed fraction of heterozygotes diagnosed with long QT syndrome. Heuristically, we found that the innate diagnostic information one learns about a variant from 3-dimensional variant location, in vitro functional data, and in silico predictors is equivalent to the diagnostic information one learns about that same variant by clinically phenotyping 10 heterozygotes. Most importantly, these data can be obtained in the absence of any clinical observations. CONCLUSIONS: We show how variant-specific features can inform a prior probability of disease for rare variants even in the absence of clinically phenotyped heterozygotes.